On 
On 



X. 



Low-dimensional dynamics embedded in a plane Poiseuille flow turbulence 
Traveling-wave solution is a saddle point? 

Sadayoshi TohQ and Tomoaki ItanoQ 
Division of Physics and Astronomy, Graduate School of Science 
Kyoto University, Kyoto 606-8502, Japan 
(May 1, 1999) 

The instability of a streak and its nonlinear evolution are investigated by direct numerical simula- 
tion (DNS) for plane Poiseuille flow at Re — 3000. It is suggested that there exists a traveling- wave 
solution (TWS). The TWS is localized around one of the two walls and notably resemble to the 
coherent structures observed in experiments and DNS so far. The phase space structure around 
this TWS is similar to a saddle point. Since the stable manifold of this TWS is extended close to 
the quasi two dimensional (Q2D) energy axis, the approaching process toward the TWS along the 
stable manifold is approximately described as the instability of the streak (Q2D flow) and the suc- 
ceeding nonlinear evolution. Bursting corresponds to the escape from the TWS along the unstable 
manifold. These manifolds constitute part of basin boundary of the turbulent state. 
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PACS numbers: 47.27.Nz,47.20.Ft,47.27.Eq,05.45.-a 



In wall-turbulent shear flows, turbulence production is 
believed to occur or at least to be triggered in the near- 
wall region in a deterministic and intermittent way, and 
also to be related to coherent structures: bursting jj]. 
Although the mechanism of the bursting has not been 
well understood, the coherence of this process suggests 
that a low-dimensional dynamics is embedded even in 
fully-developed turbulence 0. This low-dimensionality 
observed in subcritical flows S is quite different from 
weakly nonlinear stage of supercritical flows such as 
Benard convection Q, because the fully-developed tur- 
bulence coexists. That is, the low-dimensional dynamics 
should be connected to the huge-dimensional dynamics 
globally. 

This global structure seems to be hopelessly compli- 
cated. However, a minimal flow unit for maintaining 
turbulence was found by Jinemez and Moin js| in plane 
Poiseuille flow. Owing to this work we can focus on the 
elementary process of turbulence in the near-wall region 
from the perspective of dynamical systems without being 
bothered by spatial coexistence of various stages. Hamil- 
ton et al. H examined the minimal flow unit of plane 
Couette flow and found that turbulence is sustained by 
not random but quasi-cyclic process passing quiescent 
and activated periods in turn. The quiescent period is rel- 
atively laminar where coherent structures, the so-called 
streaks, are formed, develop and become unstable. Then 
their instability triggers an activated or turbulent period. 
They called this process self-sustaining process (SSP). 
Note that the Reynolds numbers (Re) used in their work 
are quite low, so they observed spatial chaos rather than 
fully-developed turbulence. 

This picture seems to work well even in realistic situ- 
ations. There are, however, still open questions. One is 
on the nonlinear evolution after the instability of streaks 
and another is why the unstable streaks are regenerated. 
The low-dimensional model introduced heuristically by 



Waleffc [Q has partially answered to these. The keys are 
the existences of steady solutions and a homoclinic orbit. 
This model, however, is wholly low-dimensional and not 
related to fully-developed turbulence. 

Recently, Waleffe obtained traveling-wave solutions 
(TWS) in plan Poiseuille flow by continuing from steady 
solutions in plane Couette flow for stress boundary con- 
ditions [0. These solutions are quite similar in shape to 
the coherent structures educed by Joeng et al. || in di- 
rect numerical simulation (DNS) , although the boundary 
conditions are different. 

DNS and experiments have revealed the existence of 
coherent structures in the near-wall region that are inti- 
mately linked to the maintenance of turbulence. In fact 
the sinuous evolution of a streak observed in DNS ap- 
pears to trigger off the formation of longitudinal coher- 
ent structures and the succeeding bursting. This process 
is also the main part of the SSP picture. Schoppa and 
Hussain |l^] showed that the instability of ejected low 
speed streaks directly generates new streamwise vortices, 
internal shear layers and arch vortices by examining sta- 
bility of a vortex-less streak |ll|. The resulting 3D vor- 
tex geometry is identical to that of the dominant coher- 
ent structures educed from the near- wall turbulence, too. 
They concluded that vortex-less streaks are the main 
breeding ground for new streamwise vortices, commonly 
accepted as dominant in turbulence production. 

In this paper, we try to elucidate the relation of the 
conceptual ingredients of the SSP picture such as the 
instability of streaks and coherent structures with TWS, 
homoclinic nature, etc. in a more realistic situation by 
means of DNS. 

The numerical scheme adopted here is basically the 
same as those of Kim et al. ]l2j. For all simulations re- 
ported, we use periodic boundary conditions in x and z 
with period lengths L x — it ~ 420 + and L z = OAir ~ 
170 + , and apply no-slip boundary conditions at y = ±1: 
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x, y and z represent the streamwise, normal to the 
walls and spanwise directions, respectively. The value 
of Reynolds number is fixed to 3000. Though somewhat 
large, this system belongs to the minimal flow unit |^|. 
The flow is maintained by a constant volume flux. As 
aliasing errors are removed by 1/2 phase-shift method, 
the effective modes are 30 x 30 Fourier modes in x and 
z, and 65 Chebychev polynomials in y. 

In Fig. |l|, the evolution of Reynolds stress integrated 
over the horizontal plane, \(u x u y )(y,t)\ is shown. Tur- 
bulent fluctuations are produced close to the walls and 
ejected into the center region. These sudden production 
and ejection of fluctuations are typical bursting processes 
p3| . Bursting on both walls seems to occur indepen- 
dently. In the lower half region, an especially quiescent 
period where the flow is laminarized and a prominent 
streak develops, lasts about ten times as long as ordi- 
nary one. Even in this period, bursting continues to oc- 
cur in the upper region. This quasi-cyclic and intermit- 
tent occurrence of bursting is reminiscent of the chaotic 
homoclinic orbit p^j . 

To describe the low-dimensional dynamics suggested, 
we introduce the 2D phase space spanned by the quasi- 
2D (Q2D) and 3D components of kinetic energy of the 
normal velocity, u y , that is, E® 2D = J v (u y i2D ) 2 dV/2 and 
£3° = J v (ul°) 2 dV/2, where V denotes the whole vol- 
ume of the computation. Here we decompose the veloc- 
ity field into two parts: Q2D flow u^ 2D (?/, z,t), and 3D 
flow u 3D (x, y, z, t). The latter is constructed by Fourier 
modes only with non-zero streamwise wave number. The 
velocity field restricted in Q2D space, is damped to the 
laminar plane Poiseuille flow monotonically for all Re 
[pTsfl . In this sense, there is no fixed point on the Q2D 
axis except for the origin. 

Hamilton et al. || defined the streak and the longitudi- 
nal vortex as the x component (u^ 20 ) of Q2D flow and its 
y, z components (u^ 2D , it^ 2D ), respectively for the mini- 
mal Couette turbulence. They also regarded 3D flow as 
turbulent components. The simplicity of these definitions 
is a little bit curious because coherent structures are ob- 
served in the near-wall region and also three-dimensional. 
These definitions, however, conceptually work well. This 
suggests that quiescent stages or generating processes of 
coherent structures are well described in a low dimen- 
sional phase space. 

To see the global behavior of the system in phase space, 
we examine the nonlinear evolutions of small 3D distur- 
bances superposed to an artificially obtained Q2D veloc- 
ity field that seems to be a prototype of a well-developed 
streak, as done by Schoppa and Hussain |l(| . They used 
an analytic streak solution without longitudinal vortex 
jn]]. In contrast, our streak solution is constructed by 
removing the x-dependent modes, i.e. the 3D compo- 
nent, from the streak-dominated velocity field obtained 
in DNS at t = 337 in Fig. |l| in order to compare the 
obtained results with the real situation. 

The initial condition is as follows: 



u(x, 0) = u^(y, z, 0) + v/i C . \ D \ l u™(x, 0), 

IK (°)H 

(0.1) 

where ||/|| = (J y pdV) 1 / 2 and u 3D (a;,0) is a solenoidal 
random vector field with a given broad spectrum and 
random phase. The following results are qualitatively 
independent of the form of « 3D (:r, 0). 

We examine the dependence of the evolutions on the 
the relative amplitude, F ac of the 3D, i.e., x-dependent 
disturbance. In Fig. |[ the evolutions of E^ 2D (t) and 
Ey B (t) for several values of F ac are shown. After initial 
transient stage where the most growing mode is selected, 
3D disturbance (Ey D (t)) seems to grow exponentially 
and then saturate to an equilibrium state in damped- 
oscillating. Finally the disturbance chooses either the 
route to laminar plane Poiseuille flow or that to turbulent 
state depending on whether the value of F ac is smaller 
or larger than the critical value. We refer to this crit- 
ical value, which is close to F* c = 9.116010224 x 10~ 8 , 
as F£ c . The passing time till the final abrupt damping 
or growth gets longer as F ac approaches F£ c . The ex- 
ponential growth in the initial stage corresponds to the 
instability of the streak, while the streak (i?Q 2D (i)) con- 
tinues to be damped. 

Figure || shows the evolutions in the phase space 
spanned by E® 2B and Ey B for several values of F ac . It 
is easy to see that a fixed point like a saddle-focus exists. 
The trajectory is oscillating but approaches monotoni- 
cally the expected fixed point. 

The difference between two solutions for different F ac 
increases exponentially with the roughly constant growth 
rate, <r+ = 0.037 even in the approaching period. This 
suggests that the evolution of i? 3D is described as a mo- 
tion around a saddle point like (dX/ dt, d(Z r + iZi) j dt) = 
(<7 + A, (— er_ + iuj)(Z r + iZi))) , where i is imaginary unit 
and cr+,cr_ are positive constants. 

The complex damping rate could be also estimated, 
but the existence of another monotonically-damping 
mode makes an estimate difficult (see Fig. ||). A rough 
estimate shows that the ratio of these two damping rates 
is about 2: the damping rates are -0.004 for the real mode 
and -0.008 for the complex one. Therefore the dimension 
of the contracted space around the fixed point is at least 
larger than 4. This dynamics is more complicated than 
Waleffe's model. 

Both the stable and unstable manifolds of the saddle 
point constitute a separatrix that separates turbulent and 
laminar states. The stable manifold is extended close to 
the Q2D axis. Because of this closeness, the stability of 
streaks, i.e., of a Q2D solution, is well recognized even in 
fully-developed wall turbulence like in SSP. 

We infer that for F£ c the trajectory reaches the fixed 
point. In physical space, the solution corresponding to 
this point has a 3D shape and moves steadily with the 
velocity v — 0.75±0.05 in the streamwise direction. Thus 
the fixed point must be TWS. In this sense, this solution 
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is not a fixed-point but a periodic orbit. Since we cannot 
yield the TWS by DNS in a strict sense anyway, we regard 
the saturated state obtained for F„ c as TWS hereafter. 

This TWS is notably resemble in shape to those ob- 
tained by Waleffe for stress boundary conditions as can 
be seen in Fig. |J, although our Re is about ten times as 
large as Waleffe's. The streamwise wave length of Wal- 
effe's solution is almost the same as ours and the spanwise 
periodicity is 1.67 times as long as ours. Note that while 
Waleffe's are symmetric with respect to the centerline of 
the channel, our TWS is confined to the lower wall (see 
Fig. |^). Our TWS is also so tall in height that it reaches 
to the upper boundary of the log-law region and is not 
localized to the near-wall region. Furthermore, in wall 
turbulence an interval between adjacent streaks is about 
100 in wall unit on average. If TWS is linked to the co- 
herent structures, there could exist TWSs with shorter 
height and/or spanwise periodicity than that of the TWS 
found by us. (See Fig. ||(c). Two ejections are observed 
on the upper wall.) This suggests that other TWSs con- 
fined to the near-wall region exist and these TWSs should 
be common in realistic turbulence. 

To compare the quiescent period of the turbulence and 
the low dimensional dynamics around the TWS, we also 
project the evolution of the turbulence onto the 2D phase 
space as shown in Fig.||. Here we only use the energies 
contained in the lower half volume because the other tur- 
bulent evolution occurs independently in the upper half 
volume (see Fig. ^|(c)). It is quite surprising that the 
projected trajectory is close to the unstable manifold, 
because wall turbulence is usually coherent or in order 
only in the near-wall region. In the bursting period the 
system goes away from the TWS along the unstable man- 
ifold. This suggests that the instability not of a streak 
but of a TWS is the origin of the bursting process. In 
Fig. H we show the solutions on the unstable manifold 
and the trajectory of turbulence marked in Fig. ^. From 
these figures, the resemblance is apparent. 

The most unstable mode of the streak has the same 
wavelength as the TWS itself for L x — it. This seems to 
imply that the instability of streaks directly breed coher- 
ent structures, i.e., TWS. This is, however, not always 
true. Indeed, when L x is doubled, the most unstable 
mode to the streak is saturated to another TWS with 
double the periodicity: the streamwise wave numbers of 
the most unstable modes for L x — ir and 2n are the same. 
Moreover, we believe that bursting is the escaping pro- 
cess from TWS along the unstable manifold. Thus the 
stability of TWS is more important than that of streak. 
The further study of the former instability should elu- 
cidate bursting, i.e., turbulence generation and also the 
determination of streak intervals in wall turbulence. 

As mentioned above, there may exist many TWSs. 
These TWSs may be connected with each others through 
turbulent or activated periods, although we have not un- 
derstood how they are. This generalized homoclinic or 
"multi-clinic" nature of wall turbulence seems to be the 
substance of coherent structures and quasi-cyclic evolu- 



tion, SSP. 

In the context of control of turbulence generation, the 
existence of separatrix supports the laminarization of tur- 
bulent flow by forcing the flow to be two-dimensional at 
least in the near-wall region e.g. by means of riblets or 
suction US 

The computations have been performed on the 
NEC/SX4 of YITP, Kyoto Univ. 
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FIG. 1. Evolution of Reynolds stress integrated over hori- 
zontal plane. Shade indicates \(u x u y )(y,t)\ > 0.5. 
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FIG. 2. Linear-log plot of evolutions of E® 2D (t) and 
El D {t) for FL = 9.116010224 x 10~ 8 (solid line), 
Fli = 9.1160104 x 10~ 8 (long-dashed line; o indicates 
t = 645), FH 1 = 9.11601 x 10~ 8 (short-dashed line), 
FlY = 9.1162 x 10~ 8 (dotted line), F^ c = 9.114 x 10~ 8 
(dash-dotted line). 
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FIG. 3. Log-log plots of evolutions of (E® 2D (t), E 3 y D {£)) in 
2D phase space. Correspondence between line and F ac is the 
same as that in Fig. 2. Thick solid line shows the evolution of 
the turbulence from t = 300(A) to t = 450(D) via t = 380(+) 
in Fig.l. o indicates t = 645 for F" . 




FIG. 4. Snapshot of TWS at t = 600 for Ff c . (a) longitu- 
dinal vorticity u> x (y,z) at z + — 398. (b) u x (x,z) at y + — 71. 
Shade indicates u x < 0.65 or lu x < —0.05. Arrows indicate 
the positions of the section. Thick solid lines are for u x — 0.4 
and 0.6. 
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FIG. 5. Mean streamwise velocity profile U(y,t). Solid 
line: TWS at t = 600 for F a J c in Fig.2. Dash-dotted line: 
laminar plane Poiseuille flow. Dashed line: a snapshot close 
to the unstable manifold at t = 650 for Fll. 
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FIG. 6. Snapshots. (a),(b): close to the unstable manifold 
(o in Fig. 3), (c),(d): in the bursting process of turbulence(+ 
in Fig. 3). (a),(c) u> x (x,y) at z + — 398. (b),(d) u x (x,z) at 
y + = 71. Shade indicates u x < 0.65 or lj x < —0.1. Thick 
solid lines are for u x — 0.4 and 0.6. The phase of (c) and (d) 
is shifted for ease of comparison. 
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